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Ch ■ Abstract 

We construct two different incompatible Poisson pencils for the Toda 
lattice by using known variables of separation proposed by Moser and by 
Sklyanin. 



1 Introduction 

A bi-Hamiltonian manifold M is a smooth manifold endowed with two compatible 
bi- vectors P , Pi such that 



>: [p ,Po] = [Po,Pi] = [Pi,Pi] = o, 

o 



where [., .] is the Schouten bracket. Such a condition assures that the linear combi- 
nation P Q — XPi is a Poisson pencil, i.e. it is a Poisson bi-vector for each A G C, and 
therefore the corresponding bracket {., .}o — A{., .}i is a pencil of Poisson brackets 

Dynamical systems on M having enough functionally independent integrals of 
motion Hi, . . . H n in involution with respect to the both Poisson brackets 

>'. {H i ,H j } = {H i ,H j } 1 = 0. (l.i; 

^ . will be called bi-integrable systems. 

The main aim of this note is to prove that any separable system is bi-integrable 
system. As an example we show that the Toda lattice is bi-integrable system with 
respect to two different Poisson pencils Pq + \P\ and Pq + XP*, which are related 
with two known families of the separated variables P,[7j. 

2 The separation of variables method 

A complete integral S(q, t) of the Hamilton- Jacobi equation 



dt \ dq 
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where q — (qi, . . . , q n ), is a solution S(q, t, ai, . . . , a n ) depending on n parameters 
a = (ai, . . . , a n ) such that 



det 



d 2 S(q,t,a) 



dqidaj 



^0. (2.2) 



For any complete integral of (12.11) solutions g« = qi(t, a, (3) and Pi = Pi(t, a, f3) of the 
Hamilton equations of motion may be found from the Jacobi equations 

ft = _ftS|^ p , = d -^l, i = l,...,n. (2.3) 

doti dqi 

According to (12.21) if we resolve a second part of the Jacobi equations with respect 
to parameters a>x, ■ ■ ■ ,a n °ne gets n independent integrals of motion 

a m = H m (p,q,t), m = l,...,n, (2.4) 

as functions on the phase space M with coordinates p, q. 

Definition 1 A dynamical system is a separable system if the corresponding com- 
plete integral S(q, t, a) has an additive form 

n 

S(q,t,a 1 ,...,a n ) = -Ht + ^ a x , ■ ■ ■ , « n )- (2-5) 

i=l 

Here the i-th component Si depends only on the i-th coordinate qi and a. 

In such a case Hamiltonian H is said to be separable and coordinates q are said 
to be separated coordinates for H, in order to stress that the possibility to find an 
additive complete integral of (12.51) depends on the choice of the coordinates. 
For the separable dynamical system we have 

d 2 d ( d \ 

S = —I —Sj = 0, for all j ^ k, (2.6) 



dq k dqj dq k \dqj 
such that the second Jacobi equations (12.31) are the separated equations 

Pj = -r^S j {qj,a u ...,a n ), or ^{pj, q j} a) = pj - — Sjfe, a) = . (2.7) 

Here j-th equation contains a pair of canonical variables pj and qj only. 

Proposition 1 For any separable dynamical system integrals of motion H m (p,q,t) 
{2-1$ are in the involution 



{H k ,H m } f = 0, k,m = 1, . . . ,n, 
with respect to the following brackets {•,•}/ on M 

{Ph Qj}f = $ij fj(P, l) i fe> Pj}f = {Qi, Qj}f = ° 5 ( 2 - 8 ) 
which depend on n arbitrary functions f\(p, q), . . . , f n {p, <?)• 
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PROOF: In fact we have to repeat the proof of the Jacobi theorem due by Liouville. 
Namely, differentiate the relations (12.41) by qj and then substitute momenta p k from 
separated equations (I2.7p to obtain 

dH m y dH m dpi _ / dH m dH m d 2 S j \ = Q 
~[ dlj V d<lj dPj dlj J 

It follows that for any H k and H m 



y dH k f 8H m + dHrn&%\ = y dH k dH m + y dH k dH m &Sj = Q 

j X 9pj V dqj d Pj dq) J 3 d Pj i)q } 3 d Pj <)/>, dq) 

Permute indexes k and m and subtract the resulting equation from the previous one 
we get 

y f rm^^dH^dHA =0 
3 V d Pj d qj d Pj dqj J 

The final assertion easily follows. 

Brackets {.,.}/ (12.81) are the Poisson brackets if and only if 

[Pf,Pf\ = 0, ( 2 -9) 

where 

^(-diag,/! /„> ^V^V ^ 



Proposition 2 If j -th function 

fiip A) = fj(Pj,Qj), j = h---,n, (2.11) 

depends only on the j-th pair of coordinates pj, qj then brackets {.,.}/ 112. 8\) are the 
Poisson brackets, which are compatible with canonical ones. 

PROOF: Substitute tensor Pf (12.101) into the equations 

[P ,Pf] = [P f ,P f ]=0, where P = 
one gets the following system of partial differential equations 

O r> U ) Jl o Jin U ! 

oq k op k dq k dp k 

for all j ^ k ^ i. The separable functions fj{pj,qj) (12. lip satisfy this system of 
equations. So, the corresponding tensor Pf (12. lOj) is the Poisson tensor, which is 
compatible with canonical tensor P . 
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Remark 1 Separated variables (p, q) are defined up to canonical transformations 



which have to preserve canonical tensor Pq and would change the second tensor Pj. 
It is clear, that freedom in the choice of the functions fj(pj,qj) is related with this 
freedom in the definition of the separated variables. 

According to the following proposition, the separation of variables method is 
closely related with the bi-Hamiltonian geometry. 

Proposition 3 Any separable dynamical system is bi-integrable system. 

PROOF: In order to get a pair of compatible Poisson brackets on the phase space 
M it is enough to postulate brackets (I2.8p 



between the known separated variables qj and pj for a given dynamical system. 
According to the Proposition 1 the corresponding integrals of motion H m (p,q,t) 
(12 .4p are in the bi-involution with respect to brackets {., .}o and {., .}/. 

In the next sections we demonstrate how this proposition works for the open and 
periodic Toda lattice. 

Remark 2 The separated variables (p, q) are canonical variables, which put the 
corresponding recursion operator iV = PjPq 1 in diagonal form. So, according to [21 
[3], a set of the separated variables (p, q) is said to be Darboux-Nijenhuis coordinates. 

Remark 3 For the stationary systems differentiate the separated equations (12 .7p 



then apply N* = P 1 Pf and substitute (I2.13P in the resulting equation one gets 



Pj = X(Pj,qj) 



{Pi-'lAf = Sijfj(pj,qj) 



{PhPj}f = {QhQj}f = °> 



(2.12) 




(2.13) 




i=i i=i 



It follows that 




j — 1, . . . ,n, 



that is, in matrix form, 



N*dH = dHF. 



(2.14) 



4 



Here n x n control matrix F with eigenvalues fx, . . . , f n is defined by 



F = J diag(/i, ...,/„) J, Jj 



OH; 



and dH is 2n x n matrix with entries di^- = dHj/dzi, where z = (q,p). 

Equation (12.141) means that the subspace spanned by covectors dHi, . . . , dH n is 
invariant with respect to iV* [3]. 

Example 1 For further use we introduce the special control matrix F, which is the 
Frobenius matrix 

/ci 1 ••• 0\ 
c 2 1 ••• 
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(2.15) 



: ••• 1 

\c„ o ■•• oy 

Here c& are coefficients of the characteristic polynomial Ajv(A) of the recursion 
operator N = PfP^ 1 : 

1 /2 ^ 

Ajv(A)= (det(JV-AI)J = A n - (dA™- 1 + . . . + c n ) = JJ(A - A,) . (2.16) 



3 Open Toda lattice 

Let us consider open Toda associated with the root system of &tf n type. The Hamilton 
function is equal to 

71-1 



where denotes the position of 2-th particle and Pi is its momenta, such that 

{Qi,Pj}o = 8ij, {Pi,Pj}o = Ui,qj}o = 0. (3.17) 
Consequently, the equations of motion read 



a- — Tl- li, — p9l 92 ^ — p9n-l 9n 



and 



Pi = e 



£ = 2, . . . , n — 1 . 



The exact solution is due to existence of a Lax matrix. Consider the L-operator 



A - Pi 



(3.18) 



and the monodromy matrix 



T(A) = L 1 (A)---L n _ 1 (A)L n (A) = ( ^ ^ ) , detT(A) = 1, (3.19) 



which depends polynomially on the parameter A 



T(A) 



A" + A 1 \ n - 1 + ... + A r 



5iA n_i + . 



B r , 



D 2 X n ~ 2 + . . . + D r 
The monodromy matrix satisfies Sklyanin's Poisson brackets: 

{T(A) f rO*)} = [r(A - /i), T(A) ® T(/i) ] , 
where r(A — fi) is the 4x4 rational r-matrix 



(3.20) 



r(A-/i) 



A — fj, 



-n, 



n 



/ 1 \ 

10 

10 

\ 1 J 



(3.21) 



Monodromy matrix T(A) is the Lax matrix for periodic Toda lattice, whereas 
the Lax matrix for open Toda lattice is equal to 



A B 




T Q (A) = KT(X) = 

The trace of the Lax matrix 

trT D (A) = A(X) = X n + H^- 1 + 



(A), 



K 



1 




Hr, 



{H i ,H j } = 0. 



(3.22) 



generates n independent integrals of motion Hi in the involution providing complete 
integrability of the system [7]. 

3.1 The Moser variables 

According to [6] we introduce the n pairs of separated variables Aj, //j, i — 1, . . . , n, 
having the standard Poisson brackets, 



{Aj, \j} — {/li, /ij} — 0, {Aj, fij} Q — 5, 



'.)■ 



(3.23) 



with the A-variables being n zeros of the polynomial A{X) and the /^-variables being 
values of the polynomial -B(A) at those zeros, 



A(Xi)=0, fi i = r ] - 1 \nB(\ 



1 1 ■ i 1 , . . . , 7i. 



(3.24) 



The Moser variables Xj may be identified with poles of the Baker- Akhiezer function 

—> 

\1/ associated with the Lax matrix T(A) f)3.19p 

T(A)* = A^, (^, a) = const, 

having some very special normalization a [8]. 
The interpolation data (13.241) plus n identities 



B(Xi)C(Xi) = detT(A) = 1 



allow us to construct the separation representation for the whole monodromy matrix 
T(A): 

A(A) = (A-A 1 )(A-A 2 )...(A-A n ), 



C(X) = -A(X)El 1 



(3.25) 



D(X) 



(A - Xi)A'(Xi) ' 
1 + B(X)C(X) 



A(X) ■ 

If we postulate the following second Poisson brackets (I2.12p 

{Xi,[ij}i = Xi5ij, {Xi,Xj}i = {fii,fij}i = 

one gets [9] 

{A{\),A(jj.)} 1 = {B(\),B(vl)} 1 = Q, 



(3.26) 



{A{X), B(fj)} 1 = -J—(jiA(\)B(ri -XA(»)B(xf 

The first bracket in (I3.26P guaranties that integrals of motion Hi f 1 3 . 2 2 j) from A(X) 
are in the bi-involution. 

Substitute polynomials A(X) and -B(A) in initial (p, g)-variables into the brackets 
(13. 26f) and solve the resulting equations to obtain the known second Poisson tensor 

m 

Th 1 r\ r\ 'JT- t— \ r~ j ^ O O 

1 ^ 6 dp i+ i dpi ^ Vl dqi dpi dq» dq{ 

t=i i=i i<j j 

The minimal characteristic polynomial of the corresponding recursion operator Nm = 
pvpenp- 1 j s e q Ua i to 

n n 

A Nm (X) = A(X) = l[(X - A,) = A" + J2 H ^' j ■ 

j=i i=i 

So, (A, /i)-coordinates are variables of separation of the action-angle type [3], i.e. the 
corresponding separated equations are trivial 

{H h Xj} = 0, i,j = 1, . . .,n. 

The Hamiltonians Hi H3.22j) from A(X) satisfy the Frobenius recursion relations 

N* M dHi = dH i+1 - Hi dH u (3.28) 

where = P ~ 1 P 1 open and H n+ i = 0, i.e. Hamiltonians Hi satisfy equation (I2.14p 
with the Frobenius matrix Ff (I2.15p . 



7 



Example 2 For the 3-particles open Toda lattice Hamiltonians Hi (I3.22p are 

Hi = -(P1+P2+P3), 

H 2 = P1P2 + PiPz + P2P3 - e qi ' q2 - e g2 ~ q3 , 

H 3 = -P1P2P3 + Pie 92 '" 3 + Pse 91 " 12 ■ 

It's obvious that H = Hf/2 — H 2 . The second Poisson tensor Pi ( 13.271) in the matrix 
form reads as 



/ 





-1 


-1 


Pi 








\ 




1 





-1 





P2 









1 


1 











P3 






-Pi 











_ e 9i-iJ2 












-P2 





e 9l-<?2 





_ e <?2-<?3 




\ 








-P3 





e 92-<?3 





/ 



(3.29) 



The control matrix F in (12.141) is the Frobenius matrix 

/ -Hi 1 \ 
FZ en = \ -H 2 1 , (3.30) 
\-H 3 / 

where coefficients q = —Hi coincide with integrals of motion. 



3.2 The Sklyanin variables 

According to [8] we can consider another set of the separated coordinates, which 
are poles of the Baker- Akhiezer function ^ associated with the Lax matrix T(A) 
(I3.19P having the standard normalization a = (0, 1). 

In this case the first half of variables is coming from (n — 1) finite roots and 
logarithm of leading coefficient of the non-diagonal entry of the monodromy matrix 

n-l 

B(\) = -e^H(\-v j ), (3.31) 

i=i 

Another half is given by 

n 

-In A(vj), j = 1, . . . , n — 1, and v n = ^jPi ■ (3.32) 



8=1 

In these separated variables other entries of T(A) read as 

tn \ n— 1 n— 1 n—l 

A + E v * IK A - u i) + E e_Uj II ^ • ( 3 - 33 ) 

and 

j=i i-^j 
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If we postulate the second Poisson brackets (12.121) 

{Vi, Uj}* = ViSij, {v h Vj}\ = {U h Uj}\ = 

one gets 

{A(X), A(»)}* = {B(X), BO*)}* = , (3.34) 

and 

{A(X),B( f i)y i = ^—(\A(\)B(n)-iiA(ti)B(\) 
X — a V 



n-l 



+ e~ u " [X + fi + Y,^ )B{X)B{fi). (3.35) 



i=l 



In initial (p, g)-variables the last bracket looks like 

{A(A),S0*)}i = ^-(Ai4(A)B(Ai)-A*A0i)S(A)) (3.36) 



n— 1 



+ e" 9 " A + /i + ^ Pi B(A)B(/i). 



8=1 

The first bracket in ( I3.34p guaranties that integrals of motion Hi f l 3 . 2 2 f) from A(X) 
are in the bi-involution. 

Substitute into the brackets (I3.34I) - (I3.36I) polynomials A(X) and B(X) in initial 
(p, g)-variables and solve the resulting equations to obtain the following Poisson 
tensor 

i Z^ e <Ls Pt g q . Qp. 2L; g q . oq . 



A 



9 9 „ „ <9 9 „ d d 
+ Pa — A — h e 91-92 — — A — h e 9 "- 2-9 " -1 — A 



<9p n <9g n dp 1 dp n dVn-\ dPn ' 

where p = Yli=i Pi * s a total momentum. 

The tensor P* is independent on q n and the minimal characteristic polynomial 
of the corresponding recursion operator N$ = P^Pq 1 is equal to 

A Ns = -e^(X + p)B(X). 

The normalized traces of the powers of N$ are integrals of motion for n — 1 particle 
open Toda lattice. As consequence the Hamiltonians Hi (13.221) from A(X) satisfy 
equation (12.141) with the following control matrix 

— ( q popen J , (3.38) 



where F^ en is the Frobenius matrix (12.151) associated with the recursion operator 
Nm for n — 1 particle open Toda lattice. 

The corresponding separated equations follow directly from the definitions of 
(u, v )-variables (13.321) 



A( Vj 



0, j = l,...,n-l 



and 



a\ = . 



(3.39) 



where A(X) = X n + Y17=i a «A n % an< ^ ai = are the values of integrals of motion. 
The first (n — 1) separated equations give rise to the equations of motion 



n-l 



{A(\), Vj } = A{ Vj ) n 



A - Vi 



n — 1, 



. Vj - Vi 



(3.40) 



which are linearized by the Abel transformation [7] 



n-l 



k=l 



-A'" 1 , 



A i_1 rfA 
^4(AT : 



j = l,...,n 



where {aj} is a basis of abelian differentials of first order on an algebraic curve 
z = A(X) corresponding to the separated equations (I3.39p . 

Remark 4 From the factorization of the monodromy matrix T(A) (13. 19ft one gets 

T n (A) = T n _i(A)L n (A), B n (X) = -e qn A n ^(X), 

where B n (X) is entry of the monodromy matrix T n (X) of n particles Toda lattice and 
A n _i(A) is entry of the monodromy matrix T n _i(A) of n — 1 particles Toda lattice. 

This implies that for the (n — 1) -particle chain the Moser variables A^ coincide 
with the n — 1 Sklyanin variables Uj, % = 1, . . . , n — 1 for the n-particle chain. 

Example 3 At n = 3 the Poisson tensor P* in the matrix form reads as 



Pi 



I 


-1 





Pi 





-Pi -p 


\ 


1 











P2 


~P2 -P 



















p 




-Pi 











_g<?l-<?2 


e 9l-<72 







~P2 





g<?l-<?2 





_g<?l-<?2 




\P1+P 


P2+P 


p 


_g<?l-<?2 


g<?l-<?2 





/ 



and at n=4 it looks like 
Pi = 






-1 


-1 





PI 








-Pl-P 


1 





-1 








P2 





-P2-P 


1 


1 














P3 


-P'i-P 























P 


-Pi 



















eii-12 





-P2 








C 11-12 





- e 92-93 










-P3 








e'J2-'33 





e 12-13 



\ 



\pi+p p 2 +p P3+P P -c9l"92 _ e 92-93+ C 9l-92 _ e 9 2 -93 
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The corresponding control matrices F in (12.141) are given by 



T-topen 



~Pl -P2-P3 

P1+P2 1 

-P1P2 + e 91 " 92 



(3.41) 



and 



j-iopen 



-P1-P2-P3-P4 0\ 

P1+P2+P3 1 

— p\pi— P1P3— P2P3+e qi ~ q2 — e q 2~ q 3 1 

pip2P3-pie 9 2-93-p 3 e'Ji-'J2 / 



(3.42) 



3.3 Periodic Toda lattice 

The trace of the Lax matrix T(A) (13.191) for periodic Toda lattice 

trT(X) = A(X) + D(X) = X n + H l X"- 1 + ---H n , {H i ,H j } = 0. (3.43) 

generates n independent integrals of motion Hi in the involution providing complete 
integrability of the periodic Toda lattice |7j. For instance the Hamilton function is 
equal to 

1 n re— 1 

h = - Pi 2 + e " +1 + e9n ~ qi ■ 

i=l i=l 

The Moser (A, /i)-variables (I3.24p are the separated variables for the open Toda 
lattice only. Nevertheless, the bi-Hamiltonian structure for the periodic Toda lattice 
may be constructed by using generic properties of the Sklyanin bracket (I3.20p [10J. 
According to [10] . in order to get the second Poisson tensor Pf er for the periodic 
Toda lattice we have to apply canonical transformation 

Pi -> Pi + e~ qi , p n -> p n + e qn . (3.44) 

to the tensor p° pen (13.271) for the open Toda lattice . In this case Hamiltonians (I3.43P 
satisfy equation (12. 14j) with the standard Frobenius matrix Ff (12.151) associated with 
the recursion operator N^ r = Pf er P _1 [TO] . 

The Sklyanin (u, w)-variables H3.31[) - (I3.32[) are the separated variables for the 
open and periodic Toda lattice simultaneously. Therefore, the one Poisson tensor 
P* determines the bi-Hamiltonian structure for the both Toda lattices, but the 
corresponding Hamiltonians satisfy equation (I2.14p with the slightly different control 
matrices Fg pen and F| er . 

Namely, for the open Toda lattice the control matrix F° pen is given by (I3.38p . For 
the periodic Toda lattice we have to change the first column of this control matrix 
F° pen (13381) by the rule: 

re-l 

FIJ = -P, i*T = (P + E P^ F T + F %Z> i = 2,...,n, (3.45) 

5=1 

where F°+" 2 = 0. 
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Remark 5 For the open and periodic Toda lattices we have two matrix equations 
for the same integrals of motion: 

NljdH = dH F M , N* s dH = dH F s 

So, we have various complimentary relations 

N* M N*dH = dH F M F s , or N s ( ) ~ 1 dHF M = dH F s . 

between integrals of motion Hi, recursion operators Ns, Nm and the corresponding 
control matrices Fm, Fs- 

Example 4 For the 3-particles periodic Toda lattice Hamiltonians Hi (13.431) are 

Hi = -(P1+P2+P3), 

H 2 = P1P2 + P1P3 + P2P3 ~ e qi ~ q2 - e q2 ~ q3 - e q3 ~ qi , 
H 3 = -P1P2P3 + Pie q2 ~ q:i + p 2 e q3 ~ qi + p 3 e qi ~ q2 . 

The second Poisson tensor (13.291) after canonical transformation (I3.44p reads as 

/ -1-1 e~ qi +p, e 93 \ 



Pi 



.per 



\ 






-1 


-1 e~ qi +p 1 





e 93 


1 





-1 e~ qi 


P2 


e <?3 


1 


1 


e~ qi 





p 3 + e q3 


-11 _ P1 


_ e -<?i 


-e~ qi 


_ e 9i-<?2 


e <?3-<?l 





-P2 


e qi ~ q2 





_ e <?2-<?3 


_ e <?3 


_ e <?3 


-p a - e q3 -e 93 " 91 








The corresponding control matrix F in (12.141) is the Frobenius matrix 



F per 
M 



Cl 


1 





c 2 





1 


c 3 









(3.46) 



where 
Cl = Pi 



-<n 



P2+P3 + e q3 , 

C2 = —P1P2 - P1P3 - P2P3 - (P2 + P3)e~ qi - (Pi + P2)e qa + e qi ~ q2 + e q2 ~ q3 - e qi ' q3 , 

c 3 = ((p 3 + e q3 )p 2 - e 92_93 )pi + (p 3 + e q3 )p 2 e' qi - p 3 e qi ~ q2 - e~ qi+q2 ~ q3 - e qi ~ q2+q3 . 

are coefficients of the minimal characteristic polynomial of the corresponding recur- 
sion operator N^ r = Pf er P . 

The same integrals of motion Hi for periodic Toda lattice satisfy equation (I2.14p 
with another recursion operator N s = P " 1 P 1 * and with the following control matrix 
d33HD-(j335]) 

/ -pi -p 2 -p 3 \ 



pper 



\ 



-Pi -P2-P3 

e ?i-?2 _ pip2 _|_ ( pi _|_ P2 )(2p 1 + 2p 2 + p 3 ) 
(e^ 2 -p 1 p 2 ){2p 1 + 2p 2 +p 3 ) 





P1+P2 1 
3 ?l-92 - pm ) 



(3.47) 
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4 Conclusion. 



For the Toda lattice associated with the root system of stf n type we present two 
different bi-Hamiltonian structures on M ~ M 2n . The introduced Poisson tensors 
Pi (l3~2TI) and (13371) are incompatible 

while the corresponding recursion operators N M = P\Pq 1 and N s = P*Pq X take 
the diagonal form in the Moser and the Sklyanin separated variables respectively. 

Associated with the tensor Pi brackets (13.261) were rewritten in the r-matrix 
form in [10J. It will be interesting to get the similar r-matrix formulation for the 
brackets (I3.34I) - (I3.36I) associated with the tensor P* . 

The research was partially supported by the RFBR grant 06-01-00140 and grant 
NSc5403.2006.1. 
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